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Short range correlations are introduced using unitary correlation method in a relativistic approach 
to the equation of state of the infinite nuclear matter in the framework of the Hartree-Fock approxi- 
mation. It is shown that the correlations give rise to an extra node in the ground-state wave- function 
in the nucleons, contrary to what happens in non-relativistic calculations with a hard core. The 
effect of the correlations in the ground state properties of the nuclear matter and neutron matter is 
studied. The nucleon effective mass and equation of state (EOS) are very sensitive to short range 
correlations. In particular, if the pion contact term is neglected a softening of the EOS is predicted. 
Correlations have also an important effect on the neutron matter EOS which presents no binding 
but only a very shallow minimum contrary to the Walecka model. 
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We discuss the role of short range correlations in a relativistic approach to the description of nuclear matter. Several 
procedures may be used to introduce short range correlations into the model wave function. All have something in 
O , common with the so called e s method or coupled-cluster expansion, introduced by Coester and Kiimmel pj and most 
D ■ elegantly and extensively developed and applied to diverse quantal systems by R. Bishop 0. 

Ph ', In this short note, we consider the unitary operator method as proposed by Villars which automatically guar- 
antees that the correlated state is normalized. The general idea of introducing short range correlations in systems 
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£S) with short range interactions exists for a long time 00 but has not been pursued for the relativistic case. 

Non-relativistic calculations based on realistic NN potentials predict equilibrium points which do not reproduce 
simultaneously the binding energy and saturation density. Either the saturation density is reproduced but the binding 
energy is too small, or the binding energy is reproduced at too high a density 0. In order to solve this problem, 
the existence of a repulsive potential or density-dependent repulsive mechanism Q is usually assumed. Due to 
Lorentz covariance and self-consistency, relativistic mean field theories include automatically contributions which 
■ are equivalent to rt-body repulsive potentials in non-relativistic approaches. The relativistic quenching of the scalar 
field provides a mechanism for saturation, though, by itself it may lead to too small an effective mass and too large 
^sO . incompressibility of nuclear matter, a situation which is encountered in the Walecka model 

In the non-relativistic case, we find a wound in the relative wave function if the vector interaction is stronger than 
the scalar interaction. However, this may not be the case if the interaction is strong enough for short distances. Then, 
the effective potential can actually turn attractive at short distances and the wave function may well have a node, as 
proposed by V. Neudatchin || and advocated by S. Moszkowski 01 j on the basis of the the quark structure of the 
nucleon. The situation with short range correlations may be more subtle than might be thought from a simple non- 
relativistic model. We show that a short range node in the relative wave function may be encountered in relativistic 
models, although it remains to be seen to what extent the relativistic description simulates the quark structure. 

In non-relativistic models the saturation arises from the interplay between a long range attraction and a short 
range repulsion, so strong that it is indispensable to take short range correlations into account. In relativistic mean 
5^ ' field models, the parameters are phenomenologically fitted to the saturation properties of nuclear matter. Although 
in this approach short range correlation effects may be accounted for, to some extent, by the model parameters, it 
is our aim to study explicitly the consequences of actual short range correlations. In a previous publication |12| we 
have discussed the effect of the correlations in the ground state properties of nuclear matter in the framework of the 
Hartree-Fock approximation using an effective Hamiltonian derived from the a — u) Walecka model. We have shown, 
for interactions mediated only by sigma and omega mesons, that the equation of state (EOS) becomes considerably 
softer when correlations are taken into account, provided the correlation function is treated variationally, always 
paying careful attention to the constraint imposed by the "healing distance" requirement. In the present note we will 
work within the same approach and will include also the exchange of pions and p-mesons. Preliminary results of the 
present work have been presented in 

We start by considering the effective Hamiltonian 

H = I Tpl(x)(-ia • V + pM) a( fll>p{g) dx + - I ^ a {x)^Jy)V a ^ lS {\x - y\)My)M*) d % dy, (1) 
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where the exchange of a, u>, p and tt mesons is taken into account, so that 

V a p„ S (r)= Y, V ^is(r) (2) 
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An interaction of the form II; 



(Si) ai g(Si) 7 5Ti • T2 



(3) 



where £,j = 0^75, describes the exchange of ir mesons. The first term in the above is the repulsive contact interaction 
and the second term is an attractive Yukawa potential. This can be rewritten in momentum space as 
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(Si) Q/ g(Ei) 7 j Ti ■ T 2 
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In equation ib is the nucleon field interacting through the scalar and vector potentials and a, (3 are the Dirac 
matrices. The equal time quantization condition for the nucleons reads, [ip a (x,t),'ipp(y,ty]+_- 

m 



the indices a and (3 refer to the spin. The field expansion for the field ip at time t=0 reads 

= 7f E [ U r$)%% + Vr{-k)cl _ % 

where U r and V r are 
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5 a p8(x — y), where 
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For free spinor fields, we have cosx(fc) = M/e(k), sinx(fc) = \k\/e(k) with e(k) = \/k 2 + M 2 . However, we will 
deal with interacting fields so that we take the ansatz cosx(fc) = M*(k)/e*(k), sinx(fc) = \k*\/e*(k), with e*(k) = 

~2 2 "* -* -* 

+ M* (fc), where fc* and M*(k) are the effective momentum and the effective mass, respectively, determined 
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self-consistently by the Hartree-Fock (HF) prescription. The vacuum | 0) is defined through c r g | 0) = - | 0) 
one-particle states are written \k,r) = c - | 0) ; two-particle and three-particle uncorrelated states are written, 



respectively as \k, r; k' , r') — -. - | 0), and \k, r; k' , r'; k", r") — c r ^ c , p c r „ p, I 0)> and so on. 

We now introduce the short range correlations through the unitary operator method. The correlated wave function 
HU is |*) = e 4f2 |$) where |$) is a Slater determinant and O is, in general, a n-body Hermitian operator, splitting 
into a 2-body part, a 3-body part, etc.. The expectation value of H is 

WHW) ($\e- iQ H e lf! |$) 
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In the present calculation, we only take into account two-body correlations. Let us denote the two-body correlated 
wave function by \k,r; k' ' r') — e in \k,r;k' , r') » fi 2 \k,r;k' ,r') where /i 2 is the short range correlation factor, the 
so-called Jastrow factor [l^J. For simplicity, we consider /12 = /(ri.2)) H.2 = T\ — r*2, and f(r) = 1 — (a + (ir) e~ JT 
where a, (3 and 7 are parameters. The choice of a real function for the unitary operator has to be supplemented by 
a normalization condition (see eq. ijSJl) which assures unitarity to leading cluster order. 

The important effect of the short range correlations is the expression for the correlated ground-state energy. Here, 
in the leading order of the cluster expansion, the interaction matrix element (k,r;k' ,r'\Vi2\k,r;k' ,r') of the HF 

expression is replaced by (k,r;k',r'\Vu + <i + h\k, r; k', r') — (k, r; k', r'\t\ + t2\k, r; k', r'), where U is the kinetic 
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energy operator of particle i. As argued by Moszkowski and Bethe ^^|, it is expected that the true ground-state 
wave function of the nucleus, containing correlations, coincides with the independent particle, or HF wave function, 
for inter particle distances r > r^, where ss 1 fm is the so-called "healing distance". This behavior is a consequence 
of the restrictions imposed by the Pauli Principle. A natural consequence of having the correlations introduced by a 
unitary operator is a normalization constraint on /(r), 



(/ 2 (r) - 1) d 3 r = 0. 
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The correlated ground state energy of symmetric nuclear matter reads 



k 2 dk [|fc|sin x(fc) + Mcosx(fc)] + 
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where A;, B^, Cj , / and J are exchange integrals defined in the appendix. In the above equation, the first term comes 
from the kinetic contribution, the second and third terms come respectively from the a and ui direct contributions to 
the correlated potential energy, the other terms arise from the exchange correlation contribution to the kinetic energy, 
and from the meson exchange contributions to the correlated potential energy. The direct term of the correlation 
contribution to the kinetic energy vanishes due to (JHJ), and ps and p s are, respectively, the baryon and the scalar 
densities. 

The couplings g ai g u , g p , g n , the meson masses, rrii, i — a, u), p, n and the three parameters specifying the short 
range correlation function, a, (3 and 7 have to be fixed. The couplings g a and g u are chosen so as to reproduce the 
ground state properties of nuclear matter. For the p and 7r-meson couplings we take the usual values g 2 /^ = 0.55 and 
f 2 /4n = 0.08. We choose m a = 550 MeV, m w = 783 MeV, m p = 770 MeV and m n = 138 MeV. The normalization 
condition JHJ determines [3. We fix a by minimizing the energy. The parameter 7 is such that it reproduces a 
reasonable healing distance, assuming that this quantity decreases as kp increases. Therefore, we assume that 7 
depends on kp according to 7 = a\ + 02 kp/kpo, where the parameters a% and 02 are conveniently chosen. 

In tables [i] and [H] we have tabulated the parameters used in our calculation together with the relative effective 
mass M*/M, the kinetic energy T/ps — M, the direct and exchange parts of the potential energy (Vd/pB and 
V e / Pb respectively) with correlation, and the correlation contribution to the kinetic energy T c /pb, all calculated at 
the saturation point. Notice that a HF calculation produces an EOS which is stiffer than the one obtained at the 
Hartree level. However, correlations reduce the effective mass and soften the EOS. In fact, the contribution of direct 
and exchange correlation terms are of the same order of magnitude of the other terms in the energy per particle. 
Moreover, the values of the couplings g a and g u which reproduce the saturation density and binding energy strongly 
depend on the correlations, being considerably reduced by their presence, which is quite a remarkable fact. Hence, 
short-range correlations cannot be disregarded. 

The correlation function f(r) is plotted in figure 1 as a function of the relative distance for two different situations: 
in one calculation the contact term coming from the pion contribution was included and the other curve was obtained 
excluding it. It is clear from both curves that the correlations give rise to an extra node in the dependence of 
the ground-state wave-function on the relative coordinate, contrary to what generally happens in non-relativistic 
calculations with a hard core, when the wave function acquires a wound. However, both curves are rather different 
behaviors: when the contact term is included the node appears very close to zero and the correlation function has a 
quite flat behavior. This is a sign that the correlation function used was not flexible enough to respond simultaneously 
to the repulsive contact interaction and to the attractive component of the interaction. 

We have computed the binding energies as function of the density for the Hartree, HF and HF+Corr and compared 
with the quark- meson-coupling model (QMC) 0|> as can be seen from fig. 2. The inclusion of correlations make 
the equation of state (EOS) softer than Hartree or HF calculations, if the contact term of the pion contribution is 
neglected. We also see that the inclusion of the p and 7r-mesons brings extra softness to the EOS. However, the curve 
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which describes the model with contact term and with the correlations taken into account shows a very stiff behavior. 
This is due to the fact that a very simplified parametrization of the correlation function was used, which was not 
flexible enough to respond to the repulsive and to attractive components of the interaction. A softer EOS around 
nuclear matter saturation density is also provided by QMC. 

In Fig. 3, we plot the effective mass versus density of nuclear matter. If correlations are included and the contact 
term is neglected the effective mass does not decrease so fast with the increase of density as in a Hartree or, even 
worse, HF calculation. This explains the softer behavior of the EOS with correlations. However, the variation of the 
effective mass with density is still smaller within the QMC model. 

Correlations also affect the behavior of the neutron matter equation of state (EoS), which is plotted in Fig. 4. 
In this calculation we include the four mesons a, u>, p and w. We perform the calculation with and without short 
range correlations and with and without the delta term of the 7r contribution. The delta term makes the equation of 
state very hard in both calculations: with and without short range correlations. We have already discussed that the 
correlation function should be more flexible to deal with the delta term. Both EoS without the delta term show an 
unrealistic behavior: either binding in the Hartree Fock calculation or a shallow minimum in the HF plus correlations 
calculation. Neutron star observation data are compatible with a zero density surface which would not be the case 
if a minimum at finite density would occur in the neutron star EoS. It is, however important to point out that the 
inclusion of correlations almost lifts this behavior of the neutron EoS typical of the Walecka model [2(| ■ 

We stress that a node occurs in the relative wave function, and the EOS becomes softer, when the energy is 
optimized with respect to variations of a. However, if a is set equal to 1, so that the node in the relative wave 
function is replaced by a wound, the EOS remains stiff. A variational treatment is therefore essential. 

We conclude that the explicit introduction of correlations has important effects. The behavior of the EOS and the 
values of the effective coupling constants are most sensitive to the presence of short range correlations, both when we 
keep and when we omit the contact term in the pion interaction. Finally, we observe that the presence of flexible short 
range correlations tends to soften the EOS. In fig. 2, the curve "HF-corr-with contact term" appears to contradict 
this statement. This is because the correlation function used was not flexible enough to respond simultaneously to 
the repulsive contact interaction and to the attractive component of the interaction. It is also clear that a richer 
parametrization of the correlation function, such as f(r) = 1 — (1 + ar + f3r 2 ) e~ ir , is required if the contact term is 
included. Work in this direction are in progress. 

It should be said that conclusions drawn from a study of this kind have only qualitative strength since the healing 
distance constraint imposed on the parameters of the correlation function is not completely unambiguous and since 
higher order terms of the cluster expansion of the expectation values have not been estimated. 
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The angular integrals are given by 
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TABLE I: Parameters of nuclear matter. We have used a density dependent parameter (HF+corr) 7 = 600 + 400 kp/kFo MeV 
for the correlation. These parameters were obtained with fixed: M = 939 MeV, m a — 550 MeV, ra w = 783 MeV, m p = 770 
MeV, = 138 MeV, g 2 p /4iv = 0.55 and fl/Air = 0.08 at k F0 = 1.3 fm" 1 with binding energy E B = e/p - M = -15.75 MeV. 
The value of 7 refers to saturation density. 





9t 


9u, 


a 


(MeV) 


7 

(MeV) 


Hartree 


11.079 


13.806 








HF 


10.432 


12.223 








HF+corr(cr + to) 


4.4559 


2.6098 


13.855 


-2252.448 


1000 


HF+corr(u + ui + p + n) 


3.1925 


2.199 


13.822 


-2258.037 


1000 


HF+corr(u + u) + p + n + 7r(contact term)) 


11.662 


12.930 


1.5435 


-497.391 


1000 



where 

F i (k > k')= I [f(r)Vi(r)f(r)\ e^'^ dr and Ci(k,k') = I (/2(r) - 1) e^'^ dr. 
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TABLE II: Ground state properties of nuclear matter at saturation density. 





M*/M 


T/p B - M 


V d / P B 


Ve/PB 


t c /pb 






(MeV) 


(MeV) 


(MeV) 


(MeV) 


Hartree 


0.540 


8.11 


-23.86 






HF 


0.515 


5.87 


-37.45 


15.83 




HF+corr (a + w) 


0.625 


15.95 


-73.12 


20.46 


19.96 


HF+corr (a + uj + p + n) 


0.645 


16.41 


-17.76 


-34.57 


20.16 


HF+corr (a + Lu + p + ir + tv (contact term)) 


0.517 


6.67 


-42.30 


11.11 


2.53 



contact term 
(MeV) 



6.21 
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FIG. 1: The correlation function f(r) for the calculation with and without the delta term of the tv contribution. 




FIG. 3: Effective mass as a function of density. 
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FIG. 4: Equation of state of neutron matter (solid line: <t + lu + p + tt + correlation (without delta), dash-dotted line: a + to + 
p + 7r + correlation (with delta), dotted line: a + ui + p + tv (without delta), dashed line: a + uj + p + tv (with delta) 



